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Nonparametric and arbitrage-free construction of call 
surfaces using /i-recovery 

P. M. Blacque-Florentin* B. Missaoui^ 


Abstract This paper is devoted to the application of an li -minimisation technique to 
construct an arbitrage-free call-option surface. We propose a nononparametric 
approach to obtaining model-free call option surfaces that are perfectly consistent with 
market quotes and free of static arbitrage. The approach is inspired from the 
compressed-sensing framework that is used in signal processing to deal with 
under-sampled signals. We address the problem of htting the call-option surface to 
sparse option data. To illustrate the methodology, we proceed to the construction of 
the whole call-price surface of the S&:P500 options, taking into account the arbitrage 
possibilities in the time direction. The resulting object is a surface free of both 
butterfly and calendar-spread arbitrage that matches the original market points. We 
then move on to an FX application, namely the HKD/USD call-option surface. 

Keywords: Option surface. Implied volatility. Arbitrage-free, Nonparametric, 

^i-minimisation. 


1 Introduction 

One of the fundamental procedures in mathematical hnance is the construction of 
the implied call surface (which is equivalent to modelling the implied volatility surface, 
by inversion of the Black-Scholes formula). In practice, such a surface is only observed 
on a hnite two-dimensional grid of maturity and strikes, i.e. at the points where market 
quotes are available. However, in order to price derivatives instruments involving values 
that fall outside the grid, it is necessary to be able to refine the resolution of the surface. 
As a hrst approximation, common practice to use spline interpolation between the points 
to recover a whole smooth surface. There are however some restrictions on the space of 
admissible solutions: call surfaces have to be subject to some no-arbitrage constraints 
across maturities and strikes (Carr-Madan|3], Davis-Hobson[7]). In that regard, spline- 
interpolation of the call surface introduces arbitrage - an undesirable feature. 
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Rather than focusing on call-surfaces, the literature has a tendency to deal with 
modelling the implied volatility surface directly. As this problem is a very important 
one, it is not surprising many approaches have been proposed in the literature to address 
the modelling of the volatility surface. In Fengler m, an arbitrage-free smoothing of the 
call surface relying on least-square minimisation with smoothness penalty is suggested. 
See also Orosi |18| for a empirical evaluation of such a methodology. Another approach 
based on penalised L^-norm was introduced in Fengler-Hin HU. In a recent article. 
Gatheral and Jacquier |12| proposed a parametrisation of the volatility surface and give 
conditions for the parametrised surface to be free of calendar-spread and butterfly-spread 
arbitrage. All these approaches allow to obtain an arbitrage-free surface, while aiming 
at being not too far from the original market quotes. 

The inherent difficulty in these approaches is that the no-arbitrage conditions on the 
implied-volatility surface are highly nonlinear, while working with call surfaces gives -as 
we shall recall later- only linear no-arbitrage constraints. 

In this article, we propose a linear-programming based methodology that recovers a 
call surface matching exactly the market quotes up to user-specified tolerance and that 
is arbitrage-free. This approach is inspired by signal processing theory; the option call 
surface at a fixed time can be seen as a two-dimensional signal sampled on an initial 
grid (the marke points), and we will be aiming at recovering it on a finer grid while still 
perfectly matching the original observations and remaining arbitrage-free. 

The classical result in sampling theory -the Nyquist-Shannon theorem- states that to 
recover a signal exactly, it is sufficient to sample at twice the highest frequency contained 
in the signal. However, in certain cases such a rate is too much for one’s purposes. In 
the signal processing community, compressed sensing has gained a lot of attention over 
the past years within several domains of application, notably MRI. A good introduction 
can be found in [B]. The idea behind compressed sensing is that one does not need to 
have that many samples provided that the signal is sparse in a given observation basis. 
In this article, we inspire ourselves from the compressed sensing framework by writing 
the call-option surface problem as a constrained /i-minimisation problem. 

Sectiondiscusses the problem implementation. In particular. Subsection 2.3 recalls 
the no-arbitrage conditions the solution has to satisfy, discusses what basis to take and 
the construction of the observation matrix. Subsection 2.1 recasts the problem in a 
nonparametric regression context, and subsection |2.2| treats of sparsity, oscillations in 
the signal and the objective function of the program. Finally, Subsection 2.4 give the 
final form of the optimisation problem and its recast as a simple linear-program. 

An example of application is provided in Section where, starting from a limited 
number of call quotes around the money of the S&:P500, we recover a finer call grid and 
the corresponding implied-volatility surface, by inverting the Black-Scholes formula. 

Applying the methodology to FX options on pegged currency (HKD/USD), we see 
that the methodology may fail to recover a sparse representation. Introducing structure¬ 
preserving functions in the observation matrix allows to circumvent this issue. This is 
the purpose of Section 
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2 Problem setting and /i-recovery 


Consider a set of maturities T. Given an underlying stock price S', a risk-free rate 
r : T rT and a dividend function q : T qt, we denote the forward and the discount 
factor respectively as 


F{T) = D{T) = 


In the rest of this chapter, in order to simplify the notations, we make the following 
assumption. 


Assumption 1. We assume below that for every maturity Ti G T, 


= kT^ 


F{Tiy 


This also entails that there is the same number of strikes to reeover on eaeh time sliee. 

tb T‘ 

As a consequence, the j strike at maturity Tj is written Kj '. When the maturity 

rp 

is clear from context, we shall just write Kj for Kj''. 


2.1 call-surface recovery as a (constrained) nonparametric regression 
problem 

From a practical point of view, perfectly matching the market mid-prices may be 
too much to ask for, or just impractical. Let us consider a few cases. First, if there 
is initial arbitrage in one or more quotes, as sometimes happens (e.g. non-convexity in 
strike), not only is there no solution to the problem to begin with, as convexity will be 
lost on a smaller grid as well, but one does not want to match the mid-price exactly. 
If one allows for a certain tolerance around the points, then a proper solution may be 
recovered. Second, from a practical point of view, it may be preferrable to accommodate 
a framework that allows a user-specified tolerance. It can be the market bid-ask spread 
(as in liquid markets), or particular trader’s needs/preferences. 

The problem we aim at solving is now a little more precise. Taking the (vectorised) 
call-surface C, we want to find a function / such that 

C{K,T) = f{K,T) + ei, 

where / is an unknown function, and ei in an error term, which we seek to control and 
bound. 

Now, while / could be virtually any (continuous) function, we decompose it into a 
basis (or concatenation of bases) (/i)^i. 


f = Y,Xifi{T,K). 
i=l 
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Remark 1. The family is also sometimes ealled a dietionary, or an aggregation 

of dietionaries, sueh as bases of polynomials, trigonometrie funetions, wavelets, ete. 

For obvious numerical reasons, it is necessary to truncate the sum: 

N 

f = Y,XiU{T,K) + e2, 

i=l 

for some N, as well as to evaluate the fi on a grid. If we write Q for the evaluation 
matrix of the fi on such a grid — whose construction we shall detail below — the initial 
relationship becomes 

C{Ki, Tj) = Qx + €1 + € 2 - 

We write e := ei + 62 the general error term, which we wish to remain within pre¬ 
specified bounds. Now is also the time to choose what decomposition basis to take. 
A natural family to be considered is polynomials in T and K, more precisely, tensor 
orthonormal polynomials in T and K. Orthonormal bases are a privileged choice for 
their numerical stability. We define two (truncated) families of polynomials in 

K and T respectively. 

1. The Pf are orthonormal on the set of all maturities 

i-Tm 

< >= / ^ Pf(t)P^(t)^r(dT) = S,J, M-) = 

2. Similarly, the Pj^ are orthonormal on the grid of all considered strikes across ma¬ 
turities. Denoting by K, the subset of strikes where no strike 

appears more than once, we define: 

< Pk,Pj, >= / Pk{k)Pl,{k)piKm = 


Let us now define the following matrices of ^ : 

( P^{Tm)Pk{K{TmA)) ••• P^{Tm)PK^{K{Tm,l)) 


Qmn — 


\Pf^{Tra)Pk{K{Tra,MK)) ••• Ptf{T^)P^^{K{T^,MK)) 


Then the matrix of tensor polynomials Q G ^il^TxMK)xiNTxNK) jg giniply defined by 
blocks as 

/ Qii • Qi Nt 


Q = 


\QMt 1 


Q Adj'Nj’, 
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2.2 Sparsity, li- and weighted-/i recoveries 

Now that we have determined what the basis and the observation matrix should be, 
we focus on the features x should guarantee, namely sparsity. We want the call-surface to 
be explained by as few parameters as possible. So we are seeking a parsimonious vector 
X G M'^, that is, a vector that has the smallest possible “/o-seminorm”: 

N 

Iloilo •= ^ ^ 

i=l 

Requiring parsimony of the solution is highly desirable, as it essentially requires that the 
solution does not overfit the data. We are aiming at obtaining a “meaningful” decompo¬ 
sition of / in the basis 

One problem remains though: problems of finding the sparsest solutions in the Iq 
sense are in general NP-hard. 

For example, NP-hardness of the /g-recovery with quadratic constraints was proved 
by Natarajan dZ]. See Elad |9l p. 13-14] for a brief explanation of the issue with Iq- 
minimisation problems. Essentially, if a problem is NP-hard, then it means that there is 
currently no efficient algorithm (i.e. an algorithm with polynomial-time complexity) to 
solve it. 

Because of that NP-hardness, a way to address it is to consider the convex relaxation 
of the problem, and try to minimise the li norm instead: 

N 

Iklli := 

i=l 

i.e. minimise the total energy of the signal. This is referred to as li-recovery. 

Remark 2. We point out that conditions guaranteeing equivalence of Iq and li solutions 
is the raison d’etre of the research area known as Compressed Sensing. There is a vast 
literature that shows that in the case of a random basis satisfying the so-called Restricted 
Isometry Property, it is possible to find equivalence of these two solutions. A rigorous 
description is beyond our scope here; we refer the interested reader to Davenport et al. 
m for an introduction, as well as a (highly non-exhaustive) list of articles of interest: 
Candes-Donoho-Sanders m Candes m, Candes-Tao m, Donoho-Elad m 

Remark 3. Choosing an li-norm means that we would rather recover a solution that 
does not allow for as many non-zero components with small amplitude as, for example, 
the l 2 -norms. However, note that, at that stage, this does not prevent us from obtaining 
a solution involving high-degree - hence highly oscillating - polynomials. The appearance 
of such high-degree polynomials can be seen as overfitting; taking financial consequences 
into account, this would entail obtaining noisy surfaces when looking at the state-price 
distribution of the underlying, or computing the local-volatility surface. 
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It turns out that weighting the li norm gives a way of enhancing sparsity and reducing 
oscillations, as mentioned in Rauhut-Ward [19]; i.e. replacing the /i-norm by a norm of 
the form 

N 

ll^ll WZl — ^ ^ Wi • Xij 
i=l 

where the Wi are predetermined weights. Rauhut-Ward consider various weights, and set¬ 
ting Wi = i, i.e. penalising polynomials of high order, provides a substantial improvement 
in reducing the oscillations in their dataset; this can be seen as a smoothing penalty in 
the solution in the context of Zi-recovery. Based on their results, we shall also considers 
weighted objective norms of the form 


N 

\\A\obj = Xi. 

i=l 

2.3 No-arbitrage constraints on the solution 

We recall here the conditions for a call-option surface to be arbitrage free and see how 
these constraints implement in our problem. For pricing and hedging purposes, obtaining 
the arbitrage-free implied-volatility surface (or equivalently the call-option surface) is a 
problem of practical importance. While no-arbitrage conditions on the implied volatility 
surface are highly nonlinear (Roper |20|1. absence of arbitrage on the call-option surface 
is guaranteed by the following set of conditions (see e.g. Fengler [10]): we denote by 
F[T) the forward value of the underlying at time T and by D(T) the value of a unit 
zero-coupon bond expiring at maturity T. 

1. Convexity: for all T, the mapping K i—?■ C{K,T) is convex. The reason is the 
following: Consider three call options C{Ki), C{K 2 ) and C{K^) such that Ki < 
K 2 < K3. One can then construct the butterfly derivative instrument 


Bf{Ki,K 2 , K 3 ) := {K 3 - K 2 )C{Ki) - {K^ - Ki)C{K 2 ) + (iF2 - K^)C{K 3 ), (1) 


whose payoff is given in Figure 
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Figure 1: The payoff of an asymmetric butterfly, with Ki = 100,iF 2 = 120 and ii's = 150. 


As one can see, this payoff is nonnegative. A breach in convexity entails that one 
can find three strikes and create such a butterfly position for a nonpositive price, 
hence an arbitrage. 


2 . Absence of Calendar-spread arbitrage: for two maturities Ti < T 2 and any strike 
denote its forward value at maturity T, namely 


KT2 


kt, im. 

F{T,y 


Then we need the mapping T 1 —>■ C{K'^, T) to be non-decreasing. To see that, take 
two strikes ,Ti) and C{K'^'^,T 2 ) with Ti < T 2 and 

^ k'^^F{T2)/F{Ti). 


Then, if one were to buy a call ,T 2 ) and sell an amount 


D{T2)F{T2)/{D{TyF{Ti)) 


of calls C{Ki,Ti). At time Ti, if the underlying value S{Ti) is less than then 
the final payoff is just T 2 ) > 0; otherwise, the position at Ti is 


C{K^\T2) 


D{T2)F{T2) 

D{Ti)F{Ti) 


{S{Ty 


FjT,) 

F{T2) 


kF), 


which is equal to a put option with maturity T 2 and strike by put-call parity. 
A put option having a nonnegative payoff, this position must also be nonnegative. 
So in either case, we need this position to be nonnegative. 


3. Monotonicity: we need 


-l/(Zl(T)F(r)) < 




(iF,T) < 0. 
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The left-hand side inequality means that an infinitesimal increase in the strike — 
which induces the same infinitesimal decrease in the final payoff — must produce 
at most this (discounted) decrease on the current price of the call. The right-hand 
side inequality means that, all else equal, a call with higher strike will always have 
a lower payoff than a call with lower strike. 

4. Payoff bounds: for S the current price of the underlying, for all K, T, we need 

{D{T)F{T)S - KD{T))+ < C{K,T) < SD{T)F{T). 

The first inequality means that the value of the call today must be nonnegative 
(since its payoff will always be nonnegative) and higher than the value of a forward 
contract with the same strike. The second inequality states that the value of a 
call with always yield less than simply holding the stock, adjusted by dividends. 
Indeed, currently having 

5. Final value: C{K, 0) = (5 — iF)"*", i.e. “the value of the is the value of its payoff”. 

2.3.1 Absence of butterfly-spread arbitrage 

We first consider the constraints that will guarantee that assumption 1 is satisfied, 
namely that the recovered finer surface is convex. In particular, at the extremal points 
of the call-surface Ki and Km^, the no-butterfly/convexity requirements become the 
following, where S denotes the price of the underlying stock: 

SD{T)F{T) - > 0 

C{Km^-i,T) - C{Km^,T) >Q. 

For iPi, this corresponds to writing the asymmetric butterfly condition where one 
call has strike zero (and so its value is just the underlying stock, adjusted by dividends; 
equivalently the discounted forward). For it corresponds to writing the asymmetric 
butterfly condition where one call has infinite strike, and so has value zero. 

As has been seen before, the (vectorised) recovered call-surface C is written 

C = Qx. 


Then, the butterfly condition becomes 


BfQx > R, 


where 


Bf = 

(Bf,i 

Bf,2 

0 \ 


V 0 

BjMt/ 
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(2) 



with 


B 




and with 


/ K2 K-, 

^ K3-K2 

0 1 

0 0 

0 0 

V 0 


0 

K2-K1 

K3-K2 
_ K4 — K2 
~ K4-K3 

1 

0 


0 


K3-K2 


Ka-K 3 

K5-K3 

K5-K4 


0 


0 

K4-K3 

K5-K4 


-1 


0 1 


1 / 


Ri = 


-SD{T)F{T) if i = 1 mod Mk 
0 otherwise. 


2.3.2 Absence of calendar-spread arbitrage 

Concerning the absence of calendar-spread arbitrage, the condition reads as follows. 
Denoting the disconnt factor np to time t by D(t) and F{t) the forward price np to time 
t, we reqnire, for all Ki, 


D{Tj)F(T,) C'(A-. ,T,)>0. 

Consider the matrix G G ^iMKxMr-iJxiMKxMT)^ Pqj. indices i G {I,-- - ,Mk}, j G 
{1, • • • , Mr - 1} and I e {!,■■■ , Mk}, 


G{{j -1) X Mk + i,l) = ^{{j-i)xMK+i=i} - ^{jxMK+i=i}- 

Then, the no arbitrage condition becomes 

GQx < 0. 

2.3.3 Absence of vertical-spread arbitrage 

Regarding the constraints —D{T) < ^ < 0, these can be rewritten as two matrix 
ineqnalities. Write H G with 

Rihj) = =3} - l{i - 1 =3}, 

and 

u,{i) = D(r,)F(r,)5i{(,_i,o mod Mk)} ' 

Then the absence of vertical spread rewrites as 

HQx < Ub. 
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2.3.4 Constraints on the bounds 

It can be checked (Fengler [TU] ) that if all other arbitrage conditions (butterfly, vertical 
spread and calendar spread) are already enforced, and under the assumption then the 
remaining restriction 

{F{Ti) - Kf)+ < C{Kj\Ti) < F{Ti) 

can simply be enforced by imposing a dimension-one constraint on the lowest point of 
the surface: 

c(KS,,ri)>(Fffi)-/c2^)+. 

As a consequence, when working on the call-surface, all no-arbitrage conditions are 
linear and can all be merged into one no-arbitrage matrix inequality 

Lx < J, (3) 

with L G ^{Mkx(4:Mt-1)+1)x{NkxNt) ^ J g ^Mkx{4Mt-1)+1^ 

We now have the elements to state the /i-recovery problem. 


2.4 Weighted Zl-recovery problem and LP equivalence 

We denote by C° the vector of original market mid-prices quotes. The only thing 
left to discussing is the form taken by the weighted li norm appearing in the objective 
function of our minimisation problem. The purpose of weighting the objective function 
is so that the optimiser will favour solutions 

Recall that by order of a polynomial, we refer to its degree plus one. A tensor 
polynomial should be penalised a function of its degree/order: the higher its degree/order, 
the more penalised it should be, and so the higher the Wi must be. For example, the 
constant polynomial has order one in strike and in maturity, hence its corresponding 
weight will be 2. The polynomial K^T has order three in strike and order 2 in maturity, 
and so its corresponding weight will be 5. Write x G "^^kxNt £qj. signal vector 
we aim at recovering and consider an element Xi thereof, with i = (y — 1) • Nk + -2, 
y G {1, • • • , A^t} and z G {1, • • • Then, the weight associated with Xi will be 

Wi = y + z. 

Based on the previous sections, the problem we are solving can be rewritten as: 


(W-Li) 


N 

min X] 

X&.'^KXNj, 

< such that 
\Ax - C°\ < e, 
^Lx < J. 


(WLl) 


There exists an easy method to solve this problem: notice it can be recast as a linear 
program by making the change of variable x = u — v, with u,v > 0. In standard form. 
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this becomes 


(LP) := 


min 

U,v£m’^K^^T 


N 


Y, Wi{ui + Vi), 
i=l 


such that 

\A{u — v) — C°\ < e, 
L{u — v) < J, 


u,v > 0, 


meaning that the recovery problem can be solved using classical linear-programming 
optimisers. We refer to Luenberg-Ye |16) for background on linear programs as well as 
related numerical algorithms. 


3 Application to the S& P500 call-surface 

Let us apply the above procedure to a small sample of mid-price call options of the 
S&PSOO priced on Oct 30. 2015, displayed in figureBefore applying our methodology, 



1.15 


Maturity 


Moneyness 


Figure 2: Original data, consisting in S&:P500 call prices at 9 strikes and 3 maturities on 
30-0ct-2015. Data obtained from the OBOE website. 


we illustrate the caveats cubic spline interpolation may produce in the next Subsection. 

3.1 Fitting splines 

Before proceeding with our methodology, we shall attempt cubic-spline interpolation. 
Cubic splines are perhaps the most straightforward way to interpolate a volatility surface, 
but they are known for oscillating and introducing arbitrage. While Figure seems to 
yield a fairly good recovery, the associated state-price density oscillates and takes negative 
values, exhibiting butterfly arbitrage. 
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Moneyness Maturity Moneyness Maturity 



Figure 3: Top-left: call-surface recovered from spline interpolation on implied volatili¬ 
ties. Top-right: Implied-volatility surface recovered from spline interpolation. Bottom: 
associated state-price density. The density takes negative values, indicating towards but¬ 
terfly arbitrage. Each slice of the interpolated surface has Mk = 104 strikes, for a total 
of Mt = 11 slices. 

3.2 Call-surface recovery 

Using the same dataset, we now proceed to using our Zi-recovery procedure. The 
code has been written in MATLAB, and we used the GUROBI HU optimiser. 

We point out that orthonormal polynomials behave badly on the boundaries of the 
interval of interest. The way we deal with it is by expanding the domain of maturities 
and strikes on which to perform the recovery, thus controlling possible explosions outside 
our original domain. No-arbitrage constraints are also enforced on the extended domain, 
so as to allow for extrapolation in the vicinity of the original domain; this last point, 
however, can be easily relaxed. 

In order to recover the call prices on the finer grid {Ki,Tj)i=i^... ^MK,j=i,- ,Mt^ applying 
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the algorithm returns the following call-surface, as well as the associated implied-volatility 
surface (Figure Q. We recover slices for 11 maturities, each consisting of 104 strikes. 
The basis has been truncated to allow tensor polynomials to have up to degree 13 in 
strike (corresponding to polynomial order in strike Nk = 14) and degree 6 in time 
(corresponding to polynomial order in strike Nt = 7). 



Moneyness Moneyness 


Figure 4: Left: call-surface recovered from weighted-^i recovery. Right: corresponding 
implied volatility surface, relative tolerance e = 5 x 10““^, Mk = 104, Nk = 14, Mt = 
11,Nt = 7. 


In particular, we check that we indeed lie withing the pre-specified bounds: 
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T=0.630137 T=0.630137 





Figure 5: Plots of the fitted Call prices (left graphs) and corresponding implied volatilities 
(right graphs) at market maturities. Green and red markers are very close and indicate 
the lower and upper tolerance at market points respectively. 

We also take a look at the sparsity of the underlying signal. For that purpose, we 
plot the signal x in the form of an Nk x Nt matrix in Figure 
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Polynomial order in maturity 

Polynomial order in strike 



Polynomial order in maturity 


Figure 6: Plot of the coefficients x obtained from the /i-recovery procedure on the S&P500 
data. On the right: the absolute values of the signal, viewed from above. 


As can be seen, the signal retains some sparsity features, where the highest coefficients 
are concentrated around polynomial indices with small degree in K and T. 

Finally, we take a look at the state-price density and local volatility surfaces. In 
our context, these quantities have a purpose for evaluating overhtting, by inspecting the 
shapes thereof. 

The local-volatility surface, in the case of neither interest rate nor dividends, is given 


by 


aioc{K, T) 



dC 

df' 


is given in Figure Figure 0 plots the quantities dC/dT, d^CjdK^ and the corre¬ 
sponding local volatility u/oc- Overhtting would typically result in a noisy local-volatility 
surface. The fact we obtain a bell-shaped curve for the state-price density, as has been 
documented in the literature (e.g. in mi suggests that the ht is reasonable. 
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Moneyness MamrUy 


Figure 7: Top-left: surface plot of the state-price density, or equivalently d‘^CjdK^. 
Top-right: surface plot of the first time derivative dCjdT. Bottom: the corresponding 
recovered local-volatility surface. 

Remark 4. The number of basis functions plays an important role in the surface recovery. 
Of course, if one truncates the basis too much, we may fail to recover a solution. If one 
does not truncate enough, one should be fine in theory, as adding functions would result 
in a solution at least as good. However, we have noted that expanding the space too much 
may lead GUROBI to a suboptimal solution, i.e. a solution whose objective function value 
is actually higher than one obtained with less basis functions available. 

4 An FX application: pegged currencies 

As described in Clark |5] , implied-volatility market quotes are usually given for strad¬ 
dles and butterflies and quoted in deltas, rather than for vanilla calls and puts and quoted 
in strikes. Conversion is directly available through Bloomberg; so our starting point shall 
be the bid/ask implied volatilities of vanilla calls given as a function of the maturity and 
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strike (rather than the Delta). 

Our currency pair of interest will be the HKD/USD, i.e. the value of Hong Kong 
Dollar in U.S. Dollars. The Hong Kong dollar has the specihcity of being pegged to the 
US Dollar, meaning that it is only allowed to move inside a hxed range from the U.S. 
Dollar. At the time of the writing, the Hong Kong dollar is such that one U.S. dollar has 
to be worth at least 7.75 HKD and at most 7.85 HKD. 

The dataset we shall consider for the rest of this section is the set of bid-ask vanilla 
implied volatilities of January 19 2016, with maturities ranging from one week to two 
years, and hve strikes per maturity, obtained from Bloomberg. From the same source, 
we also obtain the Hong Kong and United States interest rates, so that we can convert 
the implied volatilities into option prices. 

In this section we shall aim at calibrating this dataset with both a parametric ap¬ 
proach from the literature and the nonparametric approach described in the previous 
section. 

For a parametric benchmark. Surface Stochastic Volatility Inspired parametrisations 
(SSVI) models present an attractive framework, and are reported to perform well on 
equity such at the S&:P500; see Gatheral-Jacquier |12j . As a benchmark, we apply their 
SSVI square-root parametrisation, which guarantees an arbitrage-free surface up to very 
long, impractical expiries. While one can see that the ht works very well overall, one 
can observe that at higher maturities, the SVI ht lies outside the bid-ask spread at the 
25-Delta Put quote. 
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T = 0.02 


T = 0.038 


T = 0.058 


T = 0.085 



Log-Strike 




T = 1.01 



I I-1—I-1-1-1“^ 
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T = 1.44 


T = 1.91 



Figure 8: Fitting of the Square-root SSVI model to the HKD/USD mid implied- 
volatilities for maturities ranging from one week to two years, using Gatheral’s R-script 
(Data source: Bloomberg). On each picture, the blue curve is the fitted SVI smile. Red 
diamonds indicate the bid market data and blue diamonds indicate the ask marked data. 
For the last two maturities and the 25-Delta Put strike, the calibrated SVI smile lies 
outside of the bid-ask spread. 


Regarding these observations, and in the case where one would need to find a solution 
guaranteed to lie within the bid-ask spread for any market-quote, there may be a case 
for applying our nonparametric methodology presented in the last chapter. 

For more than eight maturities, our previous algorithm fails to converge. 

Even with only eight maturities, when inspecting the state-price density (Figure]^, 
it is extremely noisy. This means that the optimiser has to use high-degree polynomials 
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in order to constrain the shape of the final solution, resulting in a very noisy state-price 
density surface. In particular, it gives probability mass to strikes which lie outside the 
pegging . Essentially, the solution can be seen as emulating a piecewise linear fitting. 
Another way to state it is that we failed to recover a sparse solution. 



Moneyness 0.9 0 Maturity i.i vo8 1.06 1.04 1.02 1 0.98 0.96 0.94 0.92 0.9 

Moneyness 


Figure 9: On the left, the recovered state-price density using a basis of tensor polyno¬ 
mials. On the right: the same picture, viewed from above: the corresponding density is 
extremely noisy, with spikes across the whole surface. 


This can be interpreted in two ways. The main one is that even with a limited amount 
of maturities, it turns out that, on this particular dataset, polynomials may just not be a 
basis in which the call option surface is sparse. Visually inspecting the density plot, one 
can see that the surface is fairly steep: from a linear part in the money, it quickly breaks 
into an almost-flat part when strikes increases to out-of-the-money: a high number of 
polynomials would be required to fit the data. 

The second may be because of the difficulties the optimiser may encounter while opti¬ 
mising simultaneously in strike and maturity. Since time interpolation across strike-levels 
makes little sense in the FX market, it may be sensible to focus on strike optimisation. 

Following these observations, we shall relax the problem in both ways, by allowing 
all slices to have different projections, and also find a basis which satisfies our sparsity 
requirements. 

4.1 Problem reformulation 

While polynomials allow one to fit any continuous function (by the Weierstrass the¬ 
orem, see e.g. ED Theorem 7.26.]), a high number of them may be required in order to 
provide a good approximation. 
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However, in our problem formulation, and unlike least-squares procedures acting 
on the call-prices, we have the liberty of choosing beforehand the projection basis. In 
particular, we can construct functions that would naturally be close to the original bid 
and ask data points, while satisfying shape constraints such as convexity. 

The modus operandi is the following. For every slice, we construct an ad-hoc set of 
functions in the following way. 

First, we take a linear interpolation (and possible extrapolation) of the bid prices 
onto an equidistant grid of strikes, where the bounds are the minimal and maximal value 
of the grid we want to recover the call prices on at the end. A priori, and if there is no 
butterfly arbitrage, this function is convex. If there is arbitrage, this function is at least 
locally convex and can be adjusted. 

We would like to be able to preserve convexity as much as possible, while having 
smooth basis functions. In order to do so, we introduce the following moUifier. 

Definition 1. A moUifier on n > 1 is an infinitely differentiable function from 
M” to M such that 


1. (j) has compact support, 

2 . 


4>{x)dx = 1, 


3. In the sense of distributions, 

lim (j)°‘{x) := lim — ) = (5(x), 

o—>-0 a—>-0 Ol 


where S denotes the Dirac function at zero: 

5{x) = 


-|-oo if X = 0, 
0 otherwise. 


Example 4.1. A standard moUifier from M to M is given by: 



where K is a constant such that fig integrates to one over M. 
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Figure 10: The standard mollifier in one dimension, (j)s '■ x ^ exp(l/(a;^ — l))l| 2 ,|<;i 

Then, in the following, we shall refer to the mollification of a function as follows. 

Definition 2 (Mollification of a function). Let / : M —?• M 6e a continuous function. We 
call /" the mollification of f with parameter a the function 

/“(x) = f f{y)4>f{x-y)dy. (4) 

Jr 

Since (/)“ G it is easy to see that /" is (see e.g. Hirsch [15] 

for standard properties of mollifiers). We also refer to Ghomi [T^ for the application of 
mollification for convex-function smoothing. 

Definitions (Structure-preserving function). By structure-preserving function, we refer 
to any -function that takes as parameters the bid-ask original market data. For 

example, the mollification of the mid-price is such a function. 

Figure displays mollifications of the mid-price linear interpolation for different 
values of a. One can indeed see that for a small enough, /" converges to / (a consequence 
from converging to the Dirac function in the sense of distributions when a —>■ 0). 
Essentially, the mollification is a weighted moving-average procedure. In particular, it 
can be seen on this particular example that for a small enough, the mollification lies 
within the bounds of the bid-ask spread. One would reach the same conclusions starting 
for example from the mid-price linear interpolation instead of the bid. 
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Figure 11: Plot of the mollified mid-price for the maturity three weeks and various values 
of a. As a increases, the mollified mid-price curve increases. The lowest curve (in black) 
represents the un-mollified original market mid curve. Green and red crosses indicate 
the market bid and ask prices respectively. 


We now have all the elements to restate the amendments to the problem from the 
previous chapter; we “disentangle” the surface and allow each slice to have its own basis 
decomposition (instead of having one tensor basis decomposition for the whole surface). 
So, the problem rewrites as follows. 

Letting denote the transpose of a vector 2 , we aim at finding 


X := *(xf\x2 


I J.T2 

1 ’ 


,x 


T2 I 
Nk\ 


TMrj 




^Nk 


that solves 


(Llpx) 


' NX X Mt 

i=l 

< such that 
|Ax - C°\ < e, 

< P, 


with A, L and P to be defined below. First, we introduce the basis matrix Q, defined 
by blocks: 


Q 


/Qi 

Q2 

0 \ 

Vo 

QmtJ 


^ xiWx) X (iVfx'X Wx) 


with all Qi G ^ Every column of Qi corresponds to the evaluation of one function 

at the various pairs of strikes and maturities present in our grids. In our context, it can 
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be either the structure-preserving function described in Definition evaluated at every 
strike present on slice i for a certain mollification parameter a, or a low-order polynomial. 

As a consequence, the {Mk x i + row of Q, when multiplied by x, will return the 
value of the call with maturity Tj and strike *. 

Then A is defined as the submatrix of Q consisting of the rows where the pair (K, T) 
corresponds to a market-quoted call option. 

L is the matrix containing all no arbitrage conditions, defined in a similar way as 
The only difference is the re-definition of Q. 

4.2 Application of the new methodology to the FX dataset 

When taking all the market maturities, with mollified bid- and mid- prices in the 
matrix Q for a range of mollification parameters, one recovers the following surfaces of 
call prices and implied volatility (Figure[I^. Again, the code is written in MATLAB, and 
we have used the GUROBI m optimiser to solve the linear program. Every recovered 
slice consists of 406 strikes. 


(U 

_o 

a 



Maturity 0.9 

Moneyness 


Figure 12: The obtained call-surface on the HKD/USD from 19/01/2016 with maturities 
between one week and two years, using the structure-preserving ^-recovery procedure. 

We also provide the slice-per-slice plots of the above fit in Figure Every sub¬ 
plot displays the implied-volatility smile with respect to the log-strike (the log of the 
strike/forward quotient), so that one can have a visual comparison with the Square-root 
SSVI fitting of Figure 
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Figure 13: The implied-volatility smiles as a function of the log-strike obtained by 
structure-preserving /i-recovery. 

We also inspect the state-price density; see Figure [T^ Concerning the distribution, 
a few points can be noted. As usual, the highest mass is concentrated around the 
money. Moreover, the value of the density decays very rapidly around the money, which 
is consistent with the HKD currency being pegged. Finally, one can notice that the 
distribution is not unimodal, but rather trimodal, with one secondary mode appearing 
out of the money, and a third, less pronounced, mode appears in the money. Note that 
the modes appear to “drift away” from the moneyness as the maturities expand. These 
observations conclude our study of this application. 
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Figure 14: The recovered state-price density surface of the HKD/USD. Right: the same 
surface, plotted from above. The presence of one main mode around the money becomes 
apparent. A second mode can be seen out of the money. A third, less marked mode 
appears in the money. Both secondary modes tend to drift away from the moneyness as 
time passes, due to the difference in interest rates in the US and Hong Kong. 


5 Conclusion 

In this article, we have developed a methodology to recover sparse call-option surfaces 
that are arbitrage free, match market prices up to desired precision, by showing that this 
problem could be set as a ^i-recovery problem, which itself could be recast as a linear 
program. 

Applying this methodology to the S&P 500 allowed us to recover the call surface, 
exhibiting a sparse decomposition and a smooth associated local-volatility surface. 

When applying this methodology to the more challenging HKD/USD call-option 
dataset, the above described methodology failed and we had to take advantage of the 
problem structure so as to recover a finer option surface. The price to pay is a loss 
in the dynamics decomposition of the implied-volatility surface into a small number of 
explanatory functions, in favour of a recovery that relies more on consistent smoothing 
of the bid-ask spread. 
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